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Abstract 

After some introductory comments relating to antiferromagnetism of crys- 
talline O2, and brief remarks on the geometry of ozone, Hartree-Fock (HF) 
theory plus second-order M0ller-Plesset (MP2) corrections are used to pre- 
dict the nuclear structure of low-lying isomers of free-space O n clusters, for 
n = 6, 8, and 12. The equilibrium nuclear-nuclear potential energy is also 
discussed in relation to the number n of oxygen atoms in the cluster. 

PACS: 31.15.Ne, 36.40. Qv 

As is well known, the O2 molecule has parallel spin electrons in its ground 
state. This results in important magnetic correlations between molecular or 
cluster units in oxygen assemblies, especially at high density, giving rise to a 
variety of 'spin-controlled' solid phases of molecular oxygen under pressure 
A remarkable result of the triplet character of oxygen molecules is the 
existence of insulating crystalline O2 in an antiferromagnetic (AF) state [2]. 
This, to our knowledge, is the only AF insulator known in nature with a 
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single atomic species present, though metallic AFs exist in the forms of Cr 
and Mn. Spin ordering has relevant consequences on electronic transport 
properties in solid oxygen. At high pressures, solid oxygen changes from an 
insulating to a metallic state ji[ , and at sufficiently low temperatures it even 
undergoes a superconducting transition jij. In the so-called red- or e-phase 
(1, @], high-pressure solid oxygen is believed to consist of O n clusters in a 
crystal lattice, with n = 4 [7j or n = 8 [8]. As to free-space ozone, the 
nuclear structure takes the form of an isosceles triangle, though we do not 
know of any crystalline form, even under high pressure, for which O3 is the 
basic building block. Turning more briefly to Oa, at least two theoretical 
studies of its low-lying free-space isomers exist 

Here, we shall use quantum-chemical methods, first of all, to study the 
nuclear structures of some low-lying isomers of the free-space clusters O n , for 
n — 6, 8, and 12. The problem of determining the free-space structure and 
other dynamical properties of oxygen clusters has been seldom addressed in 



the literature, one exception being Ref. [Ill], where however the geometry was 
restricted to the ring shape. Here, Hartree-Fock plus second-order M0ller- 
Plesset (HF+MP2) theory 12] was applied first of all for the optimization 
of the clusters with subsequent additions of more refined methods. All cal- 
culations were performed by means of Gaussian 09 software 13), where the 



6-311 basis set [14J, [15[ was used, supplemented by polarization and diffuse 
functions (6-311++G**). 

Fig. [T] shows a low- lying isomer of 0$, which appears as a linkage be- 
tween two somewhat perturbed, of course, ozone molecules. Table [1] collects 
the various bond lengths and relevant angles associated with this nuclear 
structure. The ground-state energy is given in atomic units (a.u.), together 
with the dissociation energy into six O atoms. 

Next, we show in the left panel of Fig. [2] a predicted nuclear structure of 
Os, in the singlet (s) state, by optimization using HF+MP2 theory. Calcula- 
tions within density functional theory (DFT) using the functional of Perdew, 
Burke, and Ernzerhof (PBE) [l6j], which has been shown by Amovilli et al. 



17] to have possibly variational validity, lead to a related structure, shown 



in the right panel of Fig. [2] (see Tab. [I] for structural details). However (see 
further comments below), the possibility that this structure represents a tran- 
sition state cannot be ruled out, as that may require really careful treatment 
of electron correlation. 

Fig. [3] depicts what we find to be a low-lying isomer of O12. 

Fig. H] shows the calculated frequencies of the normal modes for the dif- 



2 



ferent oxygen clusters under consideration. While all frequencies are real in 
the case n = 6 (Fig. [TJ, some frequencies turn out to be imaginary in the case 
n = 8. This has been confirmed by further calculations at the CCSD level of 
approximation (employing the 6-311++** basis set). Specifically, while the 
Os cluster in the triplet configuration does not form at all, one finds a Og 
cluster in the singlet configuration, with a geometry much similar to the one 
shown in the left panel of Fig. [2] (see Tab. [TJfor structural details). However, 
two frequencies turn out to be imaginary, thus showing that this structure 
might correspond to a second order saddle point in the configurational energy 
(null gradient, but two imaginary frequencies). 

Then, as a further related property, we have also extracted from the 
present bond lengths, the bare electrostatic nuclear-nuclear potential energy, 
U nn , for all the low-lying isomers considered here. This is stimulated by 
the early work of Mucci and March 18| , who studied the equilibrium values 
of U nn for a wide variety of tetrahedral and octahedral molecules. Their 
surprising conclusion was that U nn correlated closely with the total number 
of electrons in these molecules. Prompted by this finding, and its relevance 
to density functional theory (DFT) |l9] , we have plotted in Fig. \5\ U nn vs the 
number n of oxygen atoms for the four clusters considered above. 

In summary, our main findings are for the nuclear arrays proposed in Fig- 
ures [TH3] for isomers of free-space O n clusters {n = 6,8, 12). An additional 
finding is that the equilibrium nuclear-nuclear potential energy of these pre- 
dicted structures accord with the conclusions reached by Mucci and March 



181 on tetrahedral and octahedral molecules. 
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Figure 1: Low- lying isomer of 06- See Table Q] for structural details, and Table [2] for 
ground-state and dissociation energy. 




Figure 2: Two low-lying isomers of Os, both in the singlet configuration: within HF+MP2 
theory (left), and within PBE theory (right). See TableQ]for structural details, and Table[2] 
for ground-state and dissociation energy. 
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Figure 3: Low-lying isomer of O12. This is essentially composed by two 06, weakly bound 
sub-clusters, separated by an average distance of 3.37 A. Within each sub-cluster, one 

finds alternating angles of 103.90° and 104.00° between two consecutive O O bonds. 

One also finds an overall singlet spin configuration. See Table [2] for ground-state and 
dissociation energy. 
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Table 1: Bond lengths (in A) and angles (in °) for the clusters 06 and Os considered 
in this work (see caption of Fig. [3] for details on the structure of O12). The last group 
of data refer to the structure of the unstable 0§ singlet cluster, obtained at the refined 
CCSD/6-311++** level of approximation. 
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Table 2: Ground-state and dissociation energies (in a.u.) for the 0„ clusters (n — 6, 8, 12) 
discussed in the text. S — s (S = t) refers to a singlet (triplet) spin configuration, 
respectively. Relevant data for other atomic and molecular clusters of interest (O n , with 
n = 1, 2, 3, 4) are also included, for reference. 
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Figure 4: Calculated frequencies of the normal modes of the three clusters under study. 
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Figure 5: Nuclear- nuclear repulsion energy, U nn (in hartrees), as a function of the number 
of electrons N in the low-lying isomers depicted in Figs. [2 [3J [3J as well as for n = 2, 3, 4. 
The two points for n — 8 refer to the two singlet clusters 0% shown in F ig. [H The solid 
line is a guide to the eye, prompted by the proposal of Mucci and March [18J. 
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